Quantum state engineering with flux-biased Josephson phase qubits by rapid adiabatic passages 
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In this paper, the scheme of quantum computing based on Stark-chirped rapid adiabatic passage (SCRAP) 
technique [Phys. Rev. Lett. 100, 113601 (2008)] is extensively applied to implement the quantum-state manip- 
ulations in the flux-biased Josephson phase qubits. The broken-parity symmetries of bound states in flux-biased 
Josephson junctions are utilized to conveniently generate the desirable Stark-shifts. Then, assisted by various 
transition pulses, universal quantum logic gates as well as arbitrary quantum-state preparations could be imple- 
mented. Compared with the usual vr-pulses operations widely used in the experiments, the adiabatic population 
passages proposed here are insensitive to the details of the applied pulses and thus the desirable population 
transfers could be satisfyingly implemented. The experimental feasibility of the proposal is also discussed. 

PACS number(s): 85.25.Cp, 03.67.Lx, 42.50.Dv. 



I. INTRODUCTION 

Quantum state engineering (QSE) based on superconduct- 
ing Josephson circuits (SJCs) jV^ has been stimulated by 
the encouraging prospects of quantum cornputing (QC) and 
quantum information processirig (QIP) yroll- The SJCs in- 
clude the charge- ||6[|, flux- 0, m, and phase qubits 19!- 
[llll as well as their variants (T^-H- Moreover, QSE with 
SJCs is also a crucial approach to investigate the fundamental 
quantum phenomena, such as geometric phase lUSll and non- 
locality llm UTIl . Note that most of the current experiments 
for QSE employ the technique that is sensitive to the exact 
design of the applied pulses. A typical example is that the 
so-called 7r-pulse is usually exactly designed to transfer the 
population from one quantum state to another However, these 
duration-sensitive operations might not be practically the most 
optimal approaches to implement the desirable quantum ma- 
nipulations with sufficiently-high fidelity and efficiency. 

Besides the exactly-designed pulse operations, adiabatic 
passage (AP) technique developed in atomic physics is also 
an efficient strategy for quantum-state control. This technique 
possesses certain significant advantages, such as high transfer 
efficiency, robustness to environment noises, and less limits 
on the designs of the operational pulses. It is well-known that 
the main APs include such as the stimulated Raman AP (STI- 
RAP) U8], Stark chirped rapid AP (SCRAP) iMUl], and 
piecewise AP (PAP) ||22|1 . etc.. Certainly, all of these methods 
are competent for various population transfers. But, certain 
limits still exist in these APs. For example, STIRAPs are not 
immune to the induced ac Stark shifts [23i1 . Also, PAP demon- 
strated in recent experiments still needs a series of ultrafast 
pulses to simultaneously control the switches of two pulses 
and phase holdings. This further requires the superb operating 
techniques. In contrast, SCRAP realized in the recent exper- 
iment 1I21I1 only needs relevant controls on the amplitudes of 
pulses and thus could be immune to the inhomogeneous level 
broadening. 

AP technique has also been used in SJCs for various appli- 
cations, (see, e.g., IH-ii]). For example, the STIRAP tech- 



nique was utilized to prepare Fock state 112511 of a nanome- 
chanical oscillator coupled to a Cooper-pair box, two-qubit 
and three-qubit entangled states in coupled Josephson flux cir- 
cuits were proposed luauTT . etc.. In particular, in Ref. ll28ll 
we proposed an effective approach, by using the SCRAP tech- 
nique, to implement the fundamental logic gates without pre- 
cisely designing the durations of the operations. In this pa- 
per, we further generalize such an idea to implement the QSE 
with flux-biased Josephson circuits, including the single-qubit 
phase shift and also the ctj; -rotation operation, as well as the 
two-qubit iSWAP gate. Our basic idea is to introduce a con- 
trollable flux-perturbation to chirp the transition frequency of 
a selected flux-biased qubit. Then, by properly controlling the 
amplitude ratio of perturbation to the transition driving, desir- 
able population transfers between the selected quantum states 
could be achieved with sufficiently-high efficiency. Based 
on these operations, various single- and two-qubit logic gates 
in these SJCs could be implemented without exactly design- 
ing the durations of the applied adiabatic pulses. Further- 
more, arbitrary superposition of the logic states could also 
been implemented, in principle. Given the flux-biased Joseph- 
son phase qubits and their manipulations have already been 
demonstrated, the present proposal should be experimentally 
feasible. 

The paper is organized as follows: In Sec. II we give a 
brief review of the usual 7r-pulse coherent excitation, and the 
SCRAP technique used in atomic physics is also included. 
Subsequently, we give a simple description of the flux-biased 
phase qubit, a detailed approach for population adiabatic 
transfers, as well as controllable construction of superposition 
state. In Sec. Ill, we discuss how to implement the two-qubit 
gate with capacitively-coupled flux-biased phase qubits by the 
proposed SCRAP technique. Finally, conclusions and discus- 
sions are given in Sec. IV. 



II. 



POPULATION TRANSFERS BETWEEN DRIVEN TWO 
LEVELS 
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Two-level system is a fundamental model in quantum 
physics. In QIP the basic information element (i.e., qubit) is 
encoded by a well-defined two-level system. There are many 
approaches to implement the population transfers between the 



two levels of the qubit. Roughly, these approaches can be 
classified into the nonadiabatic- and adiabatic passages llisll . 



A. Nonadiabatic population passages: 7r-pulse coherent 
excitations 

The Hamiltonian of a single qubit driven by an external 
field can be generally written as (fi, = 1) 



Hs{t)^'^cj,+R{t)a^ 



(1) 



where Gz and (Tx refer to the Pauli spin operators and cjq 
the transition frequency of the qubit. R{t) = il{t) cos{vt) 
is the driving term with v being the pulse frequency and 
il{t) = e{t)fj,/h the Rabi frequency. Here, e{t) is the am- 
plitude of the applied pump pulse, and /z the matrix element 
of electric dipole moment. Suppose that the applied pump 
Iac{t) — £{t) cos{vt) is resonant with the qubit, i.e., v = wq. 
In the interaction picture, Eq. ([TJ reduces to 



Hint{t) = exp(itwoo-z/2)i?(t)crxexp(-itwocrz/2) 



(2) 



under the usual rotating-wave approximation (RWA). Corre- 
spondingly, the solution to the Schrodinger equation 
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dt 



Hint{t)Uint{t), 
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can be expressed as 



Uint{t) = CXp ( -«/ H^rit 
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Here, A{t) = /^ VL{t')dt' and I is an unit matrix. This implies 
that, if the qubit is initially prepared in the ground state, then at 
the time t the probability for the qubit evolving to the excited 
state is Pe{t) = (1 — cosA{t))/2. Therefore, in order to 
realize complete population inversion, the pulse area must be 
precisely designed as tt. Any deviation from the precise pulse 
area may result in dynamical error for the desirable population 
inversion. 



B. Adiabatic population passages: Stark Chirps 

For loosening the above rigorous requirement on exactly 
designing pulse area and improving the operational reliability, 
we add a controllable perturbation to the Hamiltonian ([TJ, i.e.. 
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In the interaction picture, the above Hamiltonian reduces to 



(5) 



H'intit) - 2 
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Note that this Hamiltonian is as the same as that in the origi- 
nal SCRAP scheme II19I1 . wherein the Hamiltonian is derived 
in the Schrodinger picture. Also, the pump pulse Iac{t) ap- 
plied in the present scheme is required to be resonant with the 
qubit. Certainly, compared with the transition frequency ujq of 
the qubit, the controllable Stark-shift term A(i) should be suf- 
ficiently small. For convenience, we rewrite the Hamiltonian 
©as 



HL^) 






n{t) 



n{t) 



2A(t) 
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where e{t) = ^/A^(t) + n^{t), and tan2z?(t) = n{t)/A{t). 
The instantaneous eigenvalues of the above time-dependent 
Hamiltonian can be straightforwardly written as: fi±{t) = 



(A(i) ± ^A^ 
vectors. 



fl'^{t))/2 with the corresponding eigen- 



\X+{t)) =sin^(i)|0) +cosi?(i)|l), (8a) 



|A_(i)) = cosi?(i) |0) - sini?(i) |1). (8b) 

In the new Hilbert space spanned by the vectors |A_(t)) and 
|A_|_(t)), the Hamiltonian (|6]l reads (see Appendix) 



^new \^ ) 



_ ( ^i^{t) -id{t) 



Under the adiabatic approximation, 

^{A\t) + n'it)f^ 



(9) 



^W^-A(t).^"W 



dt ^ ' dt 

Hamiltonian Eq. (|9]) can be further simplified to 



Had{t) = 



M-(t) 
fi+{t) 



(10) 



(11) 



The vanished nondiagonal elements denote that there is not 
any transition between the two instantaneous eigenstates 
|A_(f)) and |A_|_(i)). This implies that the qubit would pas- 
sage individually along one of the two adiabatic paths, as long 
as ?9(f) changes slowly. As a consequence, the generic solu- 
tion of the system takes the form 

«+(0)cxp(-j/gV+(i')dt')|A+(i)). (12) 

Obviously, although the population of an adiabatic state is 
conservational for no coupling between adiabatic states, the 
components in the states |0) and |1) can still vary with the 
time-dependent ^(i). In principle, one can realize arbitrary 
population distributions in |0) and |1) along one selected adi- 
abatic path. As an obvious advantage, the population transfer 
presented here is not sensitive to the pulse area. 



When the pump pulse is absent, the Hamiltonian Eq. ^ be- 
comes iJ. = A(t)|l)(l|. This indicates that a Stark-chirping 
pulse is sufficient to produce a phase shift gate Uz{a) — 
exp(m|l)(l|) with a = - J^^ A{t')dt': 



|0> ^ |0), |1> 



1 



(13) 



This is similar to the idea by lowering the potential to imple- 
ment the fast qubit's readout |[29ll . This is because that the 
Stark pulse does not destruct the population distributions in 
the states |0) and |1), but just leads to the phase accumula- 
tions 1 30']. By combining the Rabi pulse for transferring the 
populations between the two levels and the Stark pulse for 
phase shift operation, one can implement arbitrary superpo- 
sition of the states |0) and |1), with controllable probabihties 
and relative phase. 

(i) Implement the tTa;-rotation operation. We design a pulse 
sequence shown in Fig.|2a): apply only the Stark pulse at the 
first for satisfying the initial condition t?(to) = 0; and then a 
pump pulse is applied but it switches off prior to Stark pulse 
for satisfying the condition 'd{tf) — n/2. This pulse sequence 
yields the following population inversion 

-tlinv '■ \ 1^4) 

i|vI/(to)) = |l>^i^|vI/(t^.)>= |0> 

along the adiabatic paths |A+(i)) and |A_(t)), respectively. 
If the qubit resides in |0) (or |1)) originally, it would evolve 
to the state |1) (or |0)) along the adiabatic path |A_(i)) (or 
|A_|_(i))). After eliminating the additional phase via a phase- 
shift gate operation H, (tt) described above, one can realize 
single-qubit NOT gate, i.e., Umot = Hz{'K)Hinv Certainly, 
this implemented process needs to be adiabatic, otherwise the 
state will evolve along one of the two Landau-Zener tunnel- 
ing paths, which suppresses the desirable complete population 
inversion between the two logic states. 

(ii) Implement the Hadamard gate. As another example, 
we set the pulses sequence as (see Fig.|2b)): Stark pulse pre- 
cedes pump pulse to obtain the initial condition i9(io) = and 
switches off prior to the pump pulse resulting in i9(t/) = 7r/4. 
This means that 



Rh 



m.)) = 10) ^^ m,)) = ^ 

|0> + |1> 



\^{t,)) = ll) ^±^ \^{tj)) = 



(15) 



V2 



This is the standard Hadamard gate operation. 



C. Physical implementation with single flux-biased phase qubit 

SJCs provide a favorable approach to implement QC due 
to its nonlinearity. Especially, flux-biased phase qubits are 
typically utilized to perform QC with superconducting cir- 
cuits gill. 

Typically, a flux-biased phase qubit is generated by a su- 
perconducting loop (of inductance L) biased by a magnetic 



\l' 



X\) 



V, 



X 



Cj 




>^^ 



■^ 



SQUID 



FIG. 1. Circuit schematic for a single flux-biased phase qubit Jsl], 
where the X symbol denotes the JJ. Iac{t) is the microwave pump 
pulse and I^ is the biased dc-current. 



flux $e and interrupted by a Josephson junction (JJ) (with ca- 
pacitance C and critical current Iq). The advantage of this 
structure is that the generated qubit can be isolated well from 
the strongly-dissipative bias leads pll]. Using Kirchhoff's 
law, the currents along all the branches of the circuit shown 
in Fig.[T]satisfy the equation 



Ij+Ic+lL=Iacit). 



(16) 



Using the Josephson current-phase relation /,/ = Iq sin S and 
the voltage-phase relation 6 = 27rVj/$o, the above equation 
can be straightforward rewritten as Il32ll 



at zn 



L 



dt = lacit). 



This equation can be further expressed by 1331 



(17) 



(18) 



with 



U{S) = -^hcos6-^Iac{t)S 
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-cos5\ -^I,,{t)S. (19) 



Here, S is gauge invariant phase difference (macroscopic vari- 



able) of the JJ, $e 

E 



■ MI^ is the applied magnetic flux, 
/o$o/27r the JJ coupling energy, and A = 2ttIqL/^q, 
(pb — 27r$ea;/$o. Also, I^ ~ Ici>o+Idc{t) is the biased Current 
with Icf,o being the constant part and Idc{t) the level-chirping 
part used for generating the Stark shift. 

The above potential function U{S) can be divided into two 
parts; time-independent and time-dependent ones, i.e.. 



U{6) = Uo{S) + Vit), 



(20) 



with 



Uq{o) = Ej I :7^ cosd 



2A 



V{t) 



$oM 
2nL 



$n 



111 



where 4>bo = 27r/0oM/$o- 

Obviously, Eq. (fTST l, i.e., the equation for the gauge in- 
variant phase difference 5, can be interpreted as that for the 
motion of a particle with mass m = Cj($o/(27r))^ mov- 
ing in the potential U{5). Certainly, the shape of this po- 
tential can be controlled by adjusting the biased current J^, 
which indirectly changes magnetic flux though the loop. The 
bounded particle moving in the potential would have discrete 
energy levels. It is well known that all bound states of natural 
atoms/molecules have definite parities, and therefore the so- 
called electric-dipole selection rule determines all the possible 
transitions between the selected levels. This rule forbids the 
transition between the states with the same parity. However, 
in certain artificial atoms generated by, e.g., the present SJCs, 
the bound states lose the definite parities, and thus the electric- 
dipole transitions between arbitrary two levels are possible. 
This provides a convenient way to design the requirable pulses 
for implementing the above population transfers. 

First, a proper magnetic flux is applied to let the junction 
has several bound levels in the potential. Usually, the lowest 
two levels with splitting-frequency wio — [Ei — EQ)/h are 
selected to encode a JJ phase qubit, and the third one |2) might 
be involved during the qubit operations. 

Second, in order to perform the expected SCRAP intro- 
duced above, a microwave pump pulse lacit) = £{t) cos wiot 
and a controllable Stark pulse Idc{t) are applied to couple the 
qubit states, and chirp the qubit's transition frequency, respec- 
tively. Under these drivings, the Hamiltonian of the above 
flux-biased JJ reads Il34ll 



with 



and 



Hs(t) = Ho + Vit), 



Hn = py2m + Un{6) = ^ E,\t){i\, 

1=0,1,2 



M 



(21) 



V{t) = -($o/27r) [—Idc{t) + Iac{t)\ S 
= -($o/27r 



M 

T 



hc{t) + Iac{t)\ Y. I*)01'^'J' 
i,j=0,l,2 



where 5ij = (*|<5|j). We assume that Stark shifts induced by 
the pump pulse are ignorable compared to those induced by 
the Stark pulse, and also that couplings (between the selected 
levels) induced by the Stark pulse are negligible compared to 
those induced by the pump pulse. In the interaction picture 
and under the usual rotating-wave appropriation, the above 



Hamiltonian can be rewritten as 

H^nt{t) - cxp(iti?o)l^(i) exp(-iti?o) 

KfSoi 

K(5io ^i(0 K(5i2exp(wt) 

K(52i cxp(— wi) ^2(0 

(22) 

where k = Vl{t)/2,v = uw - W21, Ai(t) = MIdc{t){5ii - 
Soo)/L, A2{t) = MIdc{t){S22 ~ Soq)/L. For clarity, we re- 
turn to the Schrodinger picture, and the above equation can be 
rewritten as 



Hit) = -^ 




kSi2 . (23) 

A2(i) + W*o/ 

With the experimental parameters Ill0ll35ll /q = 8.351 /lA, 
Cj = 1.2 pF L = 168 pH, L/M = 81, and I^o = 923.7 ^lA, 
one can numerically confirm that four bound states (levels) 
exist in the left well of the potential U{S), and also Sqq ~ 
1.571, (5ii = 1.598, ^22 = 1.633, 5oi = Sw = 0.076, S12 = 
S21 = 0.109, 5o2 = S20 = -0.006, ujiq/2tt = 10.981 GHz, 
uj2i/'2tt= 10.340 GHz. 

Fig. |2] shows the desirable single qubit operation, (c) de- 
notes that, under the pulse sequence in (a), the population in 
the state |1) can be completely inverted to the state |0). In- 
versely, if the state is initially prepared at the state |0), then 
this pulse sequence will drive the system evolving completely 
to the state |1). Clearly, during this SCRAP for implementing 
the cr^; -rotation, leakage to |2) (see the dashed red line in (c)) is 
really significantly small and thus could be neglected. These 
numerical results also confirm that the desirable population 
inversion can be finished within the time interval T > 10 
ns, which is really rapid compared to the sufficiently-long de- 
coherence time (typically, e.g., 120 ns 130(1 ). Analogously, 
Fig. [2jd) numerically confirms that the desirable Hadamard 
gate operation can be demonstrated by applying the pulse 
sequence in (b). Here, |1) is assumed to be populated ini- 
tially, then under the designed pulse sequence the population 
is passed to the final state ( |0) + 1 1) )/\/2 according to Eq. ( fTsT l. 

The SCRAP-based population transfer approach proposed 
here can be directly utilized to implement the readout of the 
qubit with significantly-high fidelity. Previously, the Joseph- 
son phase qubit is read out by applying a readout pulse to fast 
lower the barrier of the potential ll29ll . The aim of this opera- 
tion is to quickly enhance the tunneling probability of the up- 
per level 1 1) for being detected. Now, our population transfer 
approach provides another way to read out the qubit. This can 
be achieved by completely transferring the population of one 
of the logic states to the readout state \R) with significantly- 
high tunneling probability for detection. This approach is sim- 
ilar to that 131 by applying a tt -pulse to resonantly drive one 
of the logic states for evolving it to the readout state \R). The 
difference in our scheme is that the duration of the 7r-pulse 
is not required to be exactly designed, and also the measure- 
ment fidelity could be significantly high. This is because that 
the population of the selected logic state has been completely 
transferred to the readout state \R) via the proposed SCRAP. 
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FIG. 2. (Color online) Population evolutions in a single flux-biased phase qubit based on the proposed SCRAP, (left) Population inversion 
operation achieved by using the adiabatic pulses (a) with the parameters Idc(t) = 5exp(— 1^/5^) nA, $1(1) — 2.98exp(— 1^/2.5^) nA. Here, 
the qubit initially resides in the state 1 1), after the designed SCRAP, |0) is fully populated. The population in the upper level (i.e., |2)) is sig- 
nificantly small (its largest value is typically less than 1%), and thus during the SCRAP the influence of this level is negligible, (right) Prepa- 
ration of the superposition state (jO) + jl))/'\/2 from the state jl) by the SCRAP with the adiabatic pulses Idc(t) = 5exp(— (t -|- 5)^/2.5^) 
nA, n(t) = 1.495exp(-tV2.5^) nA. 



III. TWO-QUBIT GATE OPERATIONS IN COUPLED 
JOSEPHSON PHASE QUBITS BY SCRAP 



For the purpose of QIP, there must be lots of qubits cou- 
pled together to form a quantum register. Fundamentally, 
any two-qubit gate assisted with arbitrary single-qubit rota- 
tion generates an universal set to produce any quantum com- 
puting circuit. In the previous section, we have shown that 
the single-qubit (Tj; -rotation, phase-shift operation, and also 
the famous Hadamard gate can be implemented in Josephson 
phase qubit by the proposed SCRAP technique. In princi- 
ple, any single-qubit rotation can be generated by combining 
the typical single-qubit operations. Now, in this section we 
show to implement a typical two-qubit gate, i.e., iSWAP one, 
with two capacitively-coupled flux-biased Josephson phase 
qubits yoj. Such a typical two-qubit gate has already demon- 
strated by using the usual tt -pulse technique OOll . The con- 
dition for exactly designing the duration of the applied pulse 
will be relaxed in our scheme. 

Without loss of generality, we consider a superconducting 
circuit formed by two capacitively-coupled flux-biased JJs. 
Also, for simplicity we assume that two junctions are iden- 
tical and possess the same energy structures (due to they are 
biased by the identical magnetic fluxes). Here, the practically- 
existing capacitive coupling could be served as the constant 



pump. An additional weak dc-current (its amplitude is time- 
dependent) is applied to one of the junction and serves as the 
required Stark pulse for chirping the levels of the qubit (see 
Fig. [3]). The Hamiltonian of this circuit is OTIl 



Hi2{t) = 



E 

fe=l,2 



-($o/27r 

fe=l,2 



Ho,k + i2n/^afpip2/C„ 
M 



_7-(2) 



-V2. 



(24) 



Here, H^k = {2n/^o?pl/{2C.j) + £.,[{5^^^ - c^b?/{2\) - 
cos(S'-'^^] describes the uncoupled k\h qubit with a renor- 
malized junction capacitance Cj = Cj{l + C) with C, = 
Cm/{Cj + Cm)- Also, Cm IS the actual coupling capac- 
itance between the two qubits, I^J{t) is the chirping cur- 
rent applied to the second junction. Furthermore, Vi = 
(27r/(f>o)^piP2/C'm represents the interaction between two 
qubits with Cm = Cj{l + C)/C being the effective coupling 
capacitance iH, and V2 = -{<^o/2TT)MI^l\t)S^^'> /L, re- 
lated to the applied chirping current, denotes the additional 
perturbation on the second qubit. 

Suppose that the chirping current is sufficiently weak, such 
that the dynamics of each qubit is still safely limited within 
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FIG. 3. Circuit diagram of the capacitively-coupled phase 
qubits 11331 • Here, two flux-biased Josephson phase qubits are cou- 
pled by the capacitance Cm- The time-dependent dc current /^^ (f ) 
is applied to the second qubit to generate the desirable Stark shift. 



the subspace 0^ = {|Ofe), |U), |2fc)}, ^^o \^k){lk\ = 1. As a 
consequence, the circuit evolves within the total Hilbert space 
%k ~ 01 ® 02- Under the usual rotating-wave appropriation, 
Vi and V2 can be rewritten as 

Vi{t) = cxp{itHoi) cxp{itHo2)Vi exp(— ifi^oi) cxp{—itHo2 

^0 ^m '- 






/_^\^j){ji\PijPji 



and 



- \02){n\pi%^^^ expi-iiu[^^t) exp{t4ft) 

- \n)mp[]lp^^^ cxpiiLoilh) expi-iJi^t) 

- |20)(llb^VPoi^ expiiiogh) expi-iJ^h) 



y2(t) = exp{itHo2)V2exp{-itHo2) 



33 



(25) 



(26) 



-m\-ls\^) 



~ —ihp'-j. Under 



with i,j = 0, 1, 2, and p^j 

the condition that Stark shifts of the levels are relatively-weak, 
the matrix elements of momenta p^-- ~ p\-' = pij. In fact, 
for the circuit with effective coupling coefficient C = 0.0017 
considered here, our numerical calculations show that: ppg = 
0.271, Ki = 0.779,p^2 = l-219,p;,i = Ko = 6.465,p;,2 = 
p^o = 1-059,^2 =P2i= 8.761. 

In the interaction picture, one can easily check that the 
dynamics of the system exists three invariant subspaces: (i) 
3i = {|00)}; (ii) Q2 = {|01),|10)}; and (iii) ^3 = 
{|02), |11), |20)}. The first subspace Qi includes only one 
quantum state |00), the reduced Hamiltonian in this subspace 
reads Hi{t) = £;oo(t)|00)(00| with 



^-W - -2^ L ^■ 



dc (0-5^0' 



27r^ 

+ ^^) r 



1 



(1) (2) 

"Poo Pm ■ 



Certainly, if the system is initially prepared in the state |00), 
then it always populates in this state. This implies that we 



have the following evolution 

|00) -^ |00). 

While, in the second subspace 3^2 = {a' = |01), b' 
the corresponding Hamiltonian can be written as 



H2it) 



with 



fla'b' 
Ab'b'it) 



^a'b' 

^b'a' Af,'b'(i) 
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Returning to the Schrodinger picture, one can get the rele- 
vant adiabatic paths and the population evolutions (Fig. Ilja, 
c)). Of course, if the two qubits are exactly resonant, i.e., 
Aft'h'(i) = 0, Hamiltonian Eq. ( l27b will lead to a periodic 
swap of the populations between the states |01) and |10) 130(1 
with the period of /i7r/(2i7a'6')- Certainly, the efficiencies 
of these population transfers are sensitive to the evolution 
time t. In order to overcome these evolution-time sensitivi- 
ties for implementing the desired swap operation, we apply 
an additional dc-current I^^ {t) = jt with 7 = 2 nA/ns to 
chirp the levels of the second qubit. If the system is initial- 
ized to be 1 10), then after the adiabatic passage (blue line in 
Fig. Eta)), |01) can be populated. This relaxes the require- 
ment of accurately designing the interaction time between the 
qubits. Yes, such a passage needs a longer time but it still 
could be finished within coherence time. In the subspace 
53 = {a = |02),5 = |ll),c = |20)}, the Hamiltonian ex- 
pressed as 



i?o(i) nabexp{—i9t) Qac 

H3{t)= \ nbaexp{iet) E[{t) Qbcexpiiet) 

^ca flcb<SXp{ — i9t) £2(1) 



(28) 



Here 



F'lt) - r*°^^7-(2)mA(2)4_r^''^2 i „(i)„(2) 

withi e {0, l,2},j = 2-i, and 

^ab = ^ba = i-^f-^;-PwP^i2^ 

*0 L^m 



0„ 



Qr 



(^\2_}_ (1) (2) 

^n^ ' r< ^20 P02 ' 

y'o ^7n 
^cb = ^6c = (-T-) ^^-7312^10' 
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Again, the corresponding Hamiltonian in the Schrodinger pic- 
ture could be written as 

E'„(t)-h0 nab nac 

H^it) = I Qba E'lit) ^nbc I . (29) 
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FIG. 4. (Color online) Population evolutions in two capacitively-coupled qubits during the designed SCRAP: the Stark pulse is generated by 
applying the weak dc-current Idc(t) = 7i with 7 = 2 nA/ns. (a) Adiabatic passages and (c) population swap betweenjOl) and |10). (b) The 
adiabatic paths and (d) the population evolution of the states in the subspace ^3 = { 1 02) , 1 11) , 1 20) } versus passage time. Here, the system 
is initially prepared in the state |11). One can see that during the passage the population in this state is unchanged (dash-dotted green), and 
also the states j02) (solid red) and |20) (dashed blue) are kept to be unoccupied. 



The adiabatic paths and the corresponding population evolu- 
tions in this subspace are also illustrated in Fig. |4|b, d). As 
what we can see that, if the population initially resides in the 
state 1 11), then it always keeps in this state during the adi- 
abatic passage designed for implementing the inversion be- 
tween the |10) and |01). Since there is not any avoid crossing 
between the state |02) (or |20)) and the state |11), the state 
|02) (or |20)) should not be populated during the above pas- 
sage. 



IV. DISCUSSIONS AND CONCLUSIONS 

In summary, we showed that populations could be adia- 
batically transferred between the selected quantum states via 
SCRAP. Typically, these transfers are insensitive to the de- 
tails of the applied pulses area. Based on these controllable 
population transfers, fundamental quantum gates and typical 
superposition states could be deterministic ally implemented. 
Our generic proposal is demonstrated by the experimentally- 
existing flux-biased Josephson phase qubits. 

Our numerical results showed that, although the designed 
passages are required to be adiabatic, the single-qubit pop- 
ulation inversion could still be finished within the duration 
as short as that of the usual 7r-pulse. However, as shown in 
Fig. (IDi that the two-qubit logic gate operation with SCRAP 



technique takes relatively long time OOll . e.g., about 300ns 
required to finish the population inversion between |01) and 
1 10) . Thus, longer decoherence is required for the coupled su- 
perconducting circuits. We also investigated the influence of 
other levels on the designed adiabatic passages. It was shown 
that the population transfers are really limited between the se- 
lected levels and the leakages to other levels are negligible. 

Note that the applied Stark pulse for chirping the levels 
should be sufficiently weak, compared with the original bias 
for defining the energy levels. For example, for realizing 
the above two-qubit operation, the maximal change ratios of 
the relevant transition frequencies (resulting from the applied 
Stark pulse) are estimated as (wj*^ — a;io)/wio = —0.617% 
and (w^Q — wio)/wio = 0.602%, respectively. Here, Lof^ 
denote the transition frequencies between the states |0) and 
|1), due to the applying of a Stark current ( ±400 nA). Also, 
the applied Stark pulses vary the transition matrix elements 
very weak, such that the influences on the actions of the usual 
pump pulses could be negligible. As a consequence, the pump 
pulses could still resonantly interact with two selected levels. 
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APPENDIX 



Suppose that the condition 

|^(i)|«A^+(i)-M-W 
is satisfied, then M[ reduces to 



(A9) 



Schrodinger equation in the interaction picture reads 



dt 



= Humit))- 



(Al) 



In the new basis of |A_(i)) and |A+(t)), 

\^it)) ^v_{t)\X^it)) +v+{t)\X+it)). (A2) 

This induces that, 

d\^{t)) 



'^f- ^ v-{t)Hum-{t))+v+{t)Hum+{t)) 



(A3) 



Multiplying the above equation left with (A_ {t)\, we get 



-zv+{t){X.{t)\-\X+{t)) 



(A4) 



Similarly, multiplying Eq. ( |A3t left with (A_|_(t)|, we have 



^d^ = «+WM+W-*«+W(A+(i)|||A+(t)) 



^iv.it){X+it)\-\X.{t)) 



This implies that 






(A5) 



(A6) 



with 

Ml = 

^_(t) - »(A_(t)|^|A_(t)) -^{X_{t)\i\X+{t)) 

-t{X+it)\i\X.{t)) M+W - *(A+(t)|^|A+(t)) 

By using the expressions of |A+(i)) and |A_(t)) in 
Eqs. 8(a,b), we have 



So, Ml can be simplified as 



(A-(()l||A,(t)) = ^ 



<^*(<i-w>-^ 



(A+(i)|^|A+(t))=0 



(A-(t)|^|A_(0)=0 



M{ 
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Thus, the solutions of Eq. 



reads 






i>+(0) 



(AlO) 



(All) 



with 



_ /cxp(-z/>_(t')d<') 







cxp(-i/oV+(^')c^i') 



Returning to the new (adiabatic) basis, the above adiabatic 
evolution solution takes the form 

\^it)) =«_(0)exp(-z/„V-(i'K)|A-(i)> + 

«+(0)exp(-z/>+(t')dt')|A+(t)). (A12) 

Here, the additional phases in adiabatic states can be regarded 
as the dynamical phases produced during the adiabatic pas- 
sages. 

Next, we show what an adiabatic condition should be satis- 
fied by controlling Stark and pump pulses. 

First, from the definition below Eq. O, 



tan2-d{t) = 



n{t) _ 2tan^ft) 
A(t) l-tan2^(i)' 



we have 



— tan2d{t) = 2{l+ian^2{){t)) — 



dt 



dt 



^A{t)~n{t)^ 



(A14) 



and 



ddji) 
dt 



dn{t 

dt 



.Ait)-n{t)^ 



2A2(i)(l + tan2 2i?(t)) 
^A{t)-n{t)^ 



2{A^t)+n^t)) 
The adiabatic condition Eq. ( |A9t implies that 



^(t)^-A(t).^"W 



dt 



dt 



«(A2(t) + r!2(i)) 



(A15) 



3/2 



(A16) 



Therefore, smooth pulses, long interaction time, and large 
Rabi frequency and detuning are needed to satisfy the desir- 
able adiabatic condition. 
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